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I. INTRODUCTION 



Materials with a significant x^ nonlinearity exhibit solitonlike beam or pulse propaga- 
tion. These materials can be modeled by the nonlinear Schrodinger equation and its variants 
[1] which are known to be integrable and possess exact soliton solutions. 

In the case of x^ materials, it has been possible to produce solitary waves through 
•^(2) .^(2) cascac ling [2]. Moreover, several particular cases of solitary wave solutions of the 
system describing this phenomena have been obtained numerically and analytically ([3-9]). 

The solitons of x^ type are known to belong to families that include periodic solutions. 
It is essential to know if this is the case for the solitary waves of the x^ type, mainly because 
the system describing them is not integrable. We do not know what properties of integrable 
systems these solitary wave solutions possess. 

Moreover, the periodic solutions containing these solitary waves as limiting cases can be 
important in their own right, specially if they are stable. For instance they may propagate 
as a background to solitary wave signals. 

In this paper we find families of periodic solutions expressed in terms of elliptic functions. 
As special limiting cases we obtain elementary trigonometric solutions and also localized 
solitary waves (including all those found earlier). We mention that not all of the solitary 
waves are stable [10-12] and that the stability of periodic solutions needs a separate study. 

In Section 2 we present the x^'-X^ cascading equations and find their Lie point sym- 
metry group. We then reduce the system to a set of four coupled real ordinary differential 
equations for traveling wave solutions. These reduced equations are further studied in Sec- 
tion 3. We show that in general the equations do not have the Painleve property [13-14] 
so the cascading equations are not integrable. We do however show that well behaved an- 
alytical solutions exist in special cases. Explicitly, elliptic function solutions are presented 
in Section 4 together with their elementary function limits. Some further special cases are 
discussed in Section 5. Conclusions are drawn in the final Section 6. 
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II. THE CASCADING EQUATIONS AND THEIR REDUCTION 



Let us write the cascading equations in the normalized form [3] 



r 

ia lt - -a lxx + a[a 2 = 0, 

2 a 2 ( 21 ) 

ia 2t - (3a 2 - i5a 2x - —a 2xx + a x = 0, 

where ai and a 2 are the (complex) envelopes of the fundamental and second harmonic wave, 
respectively. They are functions of two variables, namely t, the normalized distance along 
the wave guide and x, the normalized transverse coordinate. The constants involved are 
r = ±1 and a, (3, S (all reals). 

To obtain explicit solutions we use the method of symmetry reduction, i.e. use Lie point 



symmetries to reduce the equations (|2 . 1| ) to ordinary differential equations. We first replace 



eq.(2.1) by four real equations, putting 

a k (x,t) = R k {x,t)e irt> ^ , = 1,2, (2.2) 

< Rk < oo, < <f) k < 2tt, 

and apply a standard algorithm to find their symmetry group [15,16,17]. In the generic case, 
when all constants in the equation are arbitrary, the symmetry group of the equations is 
just three-dimensional, consisting of translations in x and t and corelated shifts of the phase 
0i — > 0i + A, 02 — > 02 + 2A. The Lie algebra of the symmetry group has a basis consisting 
of the operators 

Pi = d x , P = d t , W = 8^+20^. (2.3) 

We shall consider solutions invariant under the subgroup generated by 

X = vP 1 + ujP + {k 1 uj-uj 1 v)W, (2.4) 

where v, to, K\ and uj\ are real constants. 

The invariant solutions then have the form 
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Ql (£) = R 1 (£) e i * 1 ®e i( - Klt - ulx \ 

(2.5) 

a 2 (f) = R 2 (Z)e i4,2 ®e 2i ( Klt -" lx \ 
£ = OJX — vt. 

The four real functions R±, R 2 , </>i, 02 satisfy the following system of coupled ordinary 
differential equations 

- ^u 2 (2R' 1 (j)' 1 + R l( f)'[) + (ruu x - v)R! 1 + R ± R 2 sin = 0, (2.6) 

- § uJ 2 (2R' 2 <p 2 + J R 2 0' 2 ') + (2acjcui - u - <Jw)i^ - R\ sin = 0, (2.7) 
- r 2 u 2 {R'{ - i?!0; 2 ) + (v - ruujRift + (feu 2 - K\)Ri 

+ R X R 2 cos = 0, (2.8) 
- f cj 2 (i? 2 ' - i?20 2 2 ) - (2acuuj 1 -v- 5uj)R 2 <p' 2 

+ (2acj 2 -2 Kl -P- 25u 1 )R 2 + i? 2 cos = 0, (2.9) 

= 2 - 201. 

It is the system of equations (2.6)-(2.9) that we wish to solve. 

III. ANALYSIS OF THE REDUCED EQUATIONS 

A. Phase locked solutions 

The system (2. 6),. ..,(2. 9) is in general quite difficult to decouple and solve. It is greatly 
simplified if we impose a supplementary restriction on the phases, namely 

= kir, k e Z, i.e. sin0 = O, cos0 = e, e 2 = 1. (3.1) 

We shall call such solutions "phase locked solutions" . Let us first assume 

rauo ^ (3.2) 

and simplify notations, putting 



rujuji — v 2auuJi — v — 8u 

A= 2 ' B= 2 ' 



2auj\ — 2ki — 25uji — (3 ruj\ — 2ki 
C = 2 , D = . 



2Ni_ „ , 27V". 
A 

The equations (2. 6),. ..,(2. 9) can be rewritten as 



M 1 = j. M 2 = -. 



(3.3) 



= (3-4) 



0' 2 = £ + ^, (3.5) 

"2 

- i?^; 2 + 2Ai2x# - £>i?! - ^i?ii? 2 = 0, (3.6) 

R' 2 ' - R 2 <\>' 2 2 + 2BR 2 <\>' 2 - CR 2 - ^R\ = 0. (3.7) 



B. Case M X M 2 + 



The phase locking condition (3.1) imposes a relation between R\ and R 2 , namely 



M 2 t 2M 1 
B + -£ = 2A+—±, (3.8) 
K 2 ti x 



so that the system (3. 4),.. .,(3. 8) is overdetermined. Expressing R±, and <p 2 in terms of 
R 2 and substituing into (3.6) and (3.8), we obtain two second order ordinary differential 
equations for R 2 . These turn out to be compatible only for Ri and R 2 constant. This case 
will be considered separately in Section 5 below. 

C. Case M 1 M 2 = 

Again R\ and R 2 are constant, unless we have 

M X = M 2 = 0. (3.9) 
Let us investigate case (3.9). We have 
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<f>' 1 = A, <f)' 2 = B = 2A. (3.10) 

We put 



A = A 2 -D, B = -^-, C = B 2 -C, D = —^, (3.11) 



and obtain a system of two ordinary differential equations 

R'( + AqR\ + B0R1R2 = 0, 



(3.12) 



i?2 + C0-R2 + DqR\ — 0. 
This system is not overdetermined since the constraint (3.8) is now simply the condition 
(3.10) on the constants, i.e. 

2av — rv — 8ruj = 0. (3.13) 
Eliminating R2 from eq.(3.12), we obtain a fourth order equation for Ri, namely 

^ +C o ^-B o D Rl + AoC = 0. (3.14) 



R\ J Ri 

If the original equations (2.1) are integrable then eq.(3.14) should have the Painleve property, 
i.e. have no movable singularities other than poles [13,14]. An algorithmic test exists [13,18] 
establishing certain properties of an equation, necessary for it to have the Painleve property. 

Thus, the general solution of eq.(3.14) must allow an expansion in the neighbourhood of 
any singular point of the form 

00 

with p a negative integer, a 7^ and three of the coefficients a k arbitrary. Then R x has a 
good chance of representing the general solution of eq.(3.14), depending on four arbitrary 
constants (one of them being £ 0) the position of the pole). The values of k for which a k are 
arbitrary (i.e. are not determined by a recursion relation), are called "resonance" values. 

Substituting the expansion (3.15) into eq.(3.14), we find p = —2, = 36/(B Dq) and 
the resonance values 
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r = -1, 6, (5±iV23)/2. (3.16) 

Thus, we have only one nonnegative integer, namely r = 6, rather than the three ones needed. 
An analysis of the obtained recursion relations shows that do, 05 are fully determined, 
a 6 is indeed free and can be chosen arbitrarily. Then a 7 and all the higher terms are fully 
determined in terms of £ an d a 6 (and of course A , B , C and D ). 

Thus, eq.(3.14) does not have the Painleve property and the cascading equations (2.1) 
are not integrable. 

The Painleve analysis does however indicates that families of "well behaved" solutions 
should exist (i.e. single valued in the neighbourhood of their movable singularities), de- 
pending on one or two free parameters rather than on four. We shall find such solutions 
below. 

An alternative procedure is to solve eq.(3.12), again under the condition (3.13), for R 2 . 
We obtain the ordinary differential equation 

2(i?2 + C0-R2X-R2" + C0-R2) — (-^2" CoR'2) 2 

+4(A + B Q R 2 )(R 2 ' + C R 2 ) 2 = 0. (3.17) 

A Painleve analysis of eq.(3.17) leads to the same conclusion as that of eq.(3.14). 

D. Introduction of the Elliptic Function Equation 

Let us look for solutions R 2 (0 satisfying eq.(3.17) and also the elliptic function equation 

R' 2 2 = 74^2 + 73^2 + 72^2 + 71^2 + 7o- (3.18) 

The compatibility of eq.(3.18) and (3.17) implies 

74 = (3.19) 

and imposes six relations between the constants in (3.18) and (3.17). These allow for the 
following solutions. 
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1. 737^0, C (C -A )^0. (3.20) 



(3.21) 



The constants 7 M in eq.(3.18) are completely specified: 

25 (Co + 2A )(Co-Aq) 
7s = g-, 72 = <^o - 2A , 71 = ^ • 

(C , o + 2A ) 2 (Co-A ) 
70 = 6Bl • 

The other amplitude, Ri(£) is given directly by the expression 

R\ = ^r (BoRl + 2(A - C )R 2 + (^-Co)(2A + C7oh (3 . 22) 
and we must require R\ to satisfy 

> o (3.23) 

in the entire range of values of R 2 , a condition to be analysed below. 

Notice that eq.(3.18) will have solutions depending on just one parameter, an integration 
constant, since the coefficients 7 M are fixed in terms of A , B , C and D . The only 
constraints on the constants in the original equations (2.1) r, a, f3 and 5 and those introduced 
in the reduction procedure (2.5) namely v, u, K\ and uj\ are (3.2), (3.13) and also (3.23). 

2. 7s ^ 0, C = 0. (3.24) 

In this case we obtain 

73 = -^, 72 = -2A , 7i = -^ (3.25) 

3 Bq 

and 70 is arbitrary. Moreover, we have 



Ri 




Thus, we obtain a two parameter family of solutions, the parameters being 7 and a constant 
arising in the integration of eq.(3.18). The constraints on the coefficients are (3.2), (3.13) 
and C = 0, i.e. 

2av 2 + 10(2^10 + (3u = (3.27) 



(we have r 2 = 1). 



3. 73^0, C = A . (3.28) 
We then have 

7s = -^p 12 = ~A , 7i = (3.29) 

and 70 arbitrary. Again, (3.18) provides a two parameter family of solutions and the con- 
straint C = A is 

3av 2 + 2ru(a - 2r) fav - + 2/V = (3.30) 

and we have 

Ri = \P^-R2. (3.31) 
V M) 

4. 73 = 0. (3.32) 
In this case eq.(3.18) reduces to 

R' 2 2 = -C R 2 + 7o . (3.33) 

The solutions are 

#1 = (3.34) 

r I -^(e-e ) 2 + /i C ^0 

i? 2 = { 4 (3.35) 
- &>) C = 0. 

E. Phase Locked Solutions for a = 

We return to eq.(2.6),...,(2.9) for a = 0, sin0 = 0, cos0 = e. From eq.(2.7) we see that 
i?2 7^ const implies 
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v + 5uj = 0. (3.36) 

Eq.(2.6) can be integrated to give 



, 2N 1 ruu!-v . , 



where N is an integration constant. 
In order to have R\ ^ we impose 

2ki + /3 + 25cji ^ (3.38) 

and obtain from eq.(2.9) that we have 

#i = [e(2«i + /9 + 25cj 1 )i? 2 ] 1/2 . (3.39) 

Finally, eq.(2.8) implies that R2 satisfies the elliptic function equation 

(Kf = -aft*- - 4A Q R1 + Stf 2 - ^y +2fa|) , (3.40) 

Since both iV and 5 are arbitrary integration constants, eq.(3.40) yields a 3 parameter family 
of solutions. 

IV. SOLUTIONS IN TERMS OF ELLIPTIC FUNCTIONS AND THEIR 

LIMITING CASES 

A. General Comments on the Elliptic Function Equation 

In Section 3 we have obtained three equations of the type 

R' 2 = faR 3 + faR 2 + faR + fa, (4.1) 
where $ are real constants and fa 7^ 0. Putting 

/ = faR (4.2) 

we obtain the equation 
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/ /2 = / 3 + /V 2 + /W + A)/t (4.3) 

We introduce the roots f 1: f 2 and / 3 of the polynomial on the right hand side of eq.(4.3) 
and rewrite this equation as 

f = U-h)U-f2)U-h)- (4.4) 

The solutions of this equation can be expressed in terms of Jacobi elliptic functions [19] if the 
three roots are distinct. The case of multiple roots leads to solutions in terms of elementary 
functions. 

Let us first consider the case of three real roots and order them to satisfy /i < / 2 < /3- 
We are only interested in real solutions. 

1. h<f<h<h 
We obtain a finite periodic solution 

/ = (/ 2 -/i)sn 2 (u,k)+fi, (4.5) 

J3- h 1 
where £ is a re& l integration constant. 

2. h<f<h = h 

A special (limiting) case of solution (4.5) is obtained for f 2 — fs, namely the solitary wave 
solution 

/ = (/ 2 -/i)tanh 2 (u)+f 1 , (4.6) 

with u as in equation (4.5). 

3. h<f2<h<f 

We obtain a singular periodic solution 

sn 2 (u, k) 
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with u and k as in eq.(4.5). 

4. /i = h < fs < f 

For fi = f 2 solution (4.7) reduces to an elementary periodic singular solution namely 

sin (u) 

5. /i < h = fs < f 

For f 2 = /3 solution (4.7) reduces to a "singular solitary wave", namely 

/ = (/ 3 -A)cotanh 2 (u)+f 1 . (4.9) 

6. /i = h = fs < f 

A triple root corresponds to a "singular algebraic solitary wave" 

The case of one real and two mutually complex conjugated roots leads to singular periodic 
solutions. 

7. /i = a + ib, f 2 = a - ib, f 3 ,a,beR, b > 
The solution is 

/3+P+ ( /3 p) (u , k) 

1 + cn(u, k) 

« = VP(e-&), A; 2 = P+ 2 a p ~ /3 , P 2 = (a-/ 3 ) 2 + & 2 . (4.12) 
In the limit 6^0 solution (4.11) reduces to the solitary wave (4.6). 

B. Explicit Solutions of the Cascading Equations for 



1. The case (3.20) 



12 



We have 

P 2 = C -2A , p 1 p 3 = -hc + 2A )(C -A ) 

2 3 (4-13) 

PoP 2 3 = ^(C + 2A )\C -A ), C ^0, C ^A 

in eq.(4.3). Using eq.(3.22) we have 

R2 = ~f, R\ = 4^[/ 2 " |(4> - Co)/ + l(A - C )(2A + Co)]. (4.14) 

A globally defined solution exists if R\ is positive for all values of /. We note that the roots 
of the polynomial defining R\ are 

f a,b = l[2(A - Co) ± y/6C (C -A )]. (4.15) 

Hence, for Co (Co — A ) < the function R\ is sign-definite and can be chosen to be positive 
definite. For C (C — A ) > a more careful analysis is required. Thus solution (4.5) 
provides a global solution if the roots satisfy one of the following relations: 



f A <fl<f<f2<fB, f A <f B <fl<f< /2, 

h<f<h<f A < f B . 

Similarly, case (4.7) and (4.11) both require 



(4.16) 



f A <f B <h< f- (4.17) 

Multiple roots occur in the two excluded cases C = and C = A Q , but also in the 
allowed case 

Co = -2A (4.18) 

For A < we obtain from eq.(4.6) 

6 1 



Ro 



Ri = 



Bo cosh 2 (v^(e-eo)) 

(4.19) 

6 14| sinh(v^4^-£o)) 



V=BoEbcosh 2 ( v /rae-G>)) 
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so that we must require B D < 0. 

Similarly, for A < eq.(4.9) provides the solution 

6 1 



R2 



B sinh 2 (v^4^-6)) 



R 6 \Ap\ cosh(yra(e-eo)) 
vAB^ si nh 2 (VW(e-eo)) 

with the requirement B D Q > 0. 

For A > it is solution (4.8) that is relevant and yields: 

6A 1 



B sin 2 (v^(f -£<,)) 



64 cos(y^(e-6)) Rn>n 
V#oA) sm 2 (v^(^-eo)) 



2. The case (3.25) : C = 

Eq.(4.3) is 

f 2 = f-2A f + 4 -Alf + l0 4 -f 

with 7 arbitrary and 

i? 2 = --|-/, R 1 = - 7 ±=(-3f + 2A ), B D 
Id§ 2y/B D 

A multiple root, namely a triple one, occurs in one case only 

J = (/ - g4) , 70 = -g^a 



and we have 



1 6 

#o (£-£o) 2 



*= - 6 1 



v/^^-eo) 2 ' 
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In all other cases two of the roots are complex and / is given in eq.(4.11),(4.12). 



3. The case (3.28) : C = A 



Eq.(4.3) in this case is 



4B 2 



(4.26) 



where 70 is arbitrary and we have 
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R2 



/, Ri = 



2B J1 1 2^D~ 
The discriminant of the cubic equation f' 2 = is 

1 6 

D = ^-(Al-3B* l0 )B 2 l0 . 



f, B D >0. 



(4.27) 



(4.28) 



For D < two roots are complex, for D > they are all real and distinct, for D = we 
have a double or triple root. The absence of a first degree term in (4.26) implies that for /j 
real the three signs of /j cannot be all the same. Moreover, if a root is equal to /, = 0, it 
must be double. 

The elementary solutions in this case are: 

A 3 

1- To — 3^7- The finite solitary wave, the singular solitary wave, the singular periodic 



solution and the singular algebraic solutions in this case are: 



R9 



D °- A ° [Stanh 2 '^ 



Bo 1 2B 



(£-&))-!], A >0 



(4.29) 



Ri — \ tt-Ri 
V B 



A[ 3co tanh 2 (^(e-eo))-l], A > 



(4.30) 



D _ POO _ \M 

V Bq 2B 



sin 



Ao < 



(4.31) 



#2 




#1 



6 



Bo(Z-Zo) 



A = 



(4.32) 



respectively. 
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2. 70 = 0. The four different real solutions in this case are: 



R2 = 2^ |A)| sech2 (^^(^ " A < (4.33) 

R2 = ~\Ao\ , ^o<0 (4.34) 

2Bq sinh 2 (V^(e-eo)) 

O A 1 

R 2 = -2^. * , A >0 (4.35) 

2fi 0s in 2 (^(e-eo)) 

and for A = we reobtain solution (4.31). 

C. Solutions for a = 

Eq.(4.3) (obtained from eq.(3.40)) in this case is 

f = f- 4A f 2 - 2B Sf - AM 2 Bl (4.36) 

where S and M are arbitrary. All possibilities for the three roots f\ occur, but we have the 
restriction 

/1/2/3 > 0. (4.37) 

We have 



Rl = _ J_/, Rl = _ e (^ + /? + 2fo,) 

For i?! to be globally defined we need / to be sign definite. This imposes the following 
restrictions on the roots for each of the solutions of section (4.1). 
Solutions (4.5) and (4.6): 

h < f < h < < h or < h < f < j 2 < U- (4-39) 
Solutions (4.7), (4.8) and (4.9): 

< fi < h < /s < f or h < f2 < < h < f- (4-40) 

Solution (4.10): 
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Solution (4.11): 



o < h = h = h < /■ 
o < h < f. 



(4.41) 



(4.42) 



V. OTHER EXPLICIT SOLUTIONS 

In Section 4 we presented solutions of eq.(2.1) satisfying a±a 2 ^ and such that |ai| and 
| a.2 1 are not constant. Let us now discuss these previously rejected solutions. 

1. Solutions with a\ = 

For a^Owe put 

5 2 5 

a 2 (x, t) = u(x, t) expi[ — {(3 H )t x\. (5.1) 

2a a 

The eq.(2.1) for a 2 (x, t), a± = reduces to the linear Schrodinger equation 

a 

iu t - —Uxx = 0. (5.2) 
For a = the solution of eq.(2.1) is 

a 2 (x,t) = e' ipt uj(0, £ = x + 5t (5.3) 
where is an arbitrary (complex) function. 

2. Solutions with Ri and R 2 constant 

We require that 4>i an( i 02 should not be constant, otherwise we obtain a\ = a 2 = 0. 
Eq.(2.6),. . . ,(2.9) for Ri and R 2 constant imply that we must have sin0 = 0, i.e. the 
solutions will be phase locked. 

Explicitly, for a ^ we have 
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2 

R 2 = -^-{ 1 l-2A ll + D) 



R\ = ^(4 7l 2 - 25 7 i + C)( 7l 2 - 2A 7l + D) 



(5.4) 

0i = 7i£ + 72 



2 = 2 7l £ + 2 72 + A;tt, e=(-l) fc 
where 7i and 72 are arbitrary constants. 

For a = we have 0i and 2 as in (5.4) 

ru 2 2 rcj? — 2n\. 
R2 = 7i + {v - rwwi)7i + J 



(5.5) 

i? 2 = -e[(u + 5w)7i - (2«i + /3 + ^1)]^ 



VI. CONCLUSION 

The symmetries used in this article are only those that exist in the generic case of the 
cascading equations (2.1), i.e. for all values of the constants r, a, (3 and 5. This "generic" 
symmetry algebra is summed up in eq.(2.3), representing space and time translations and a 
shift in the phases of the fundamental and second harmonic waves. 

The reduction to the system of ordinary differential equations (2. 6). ..(2. 9) was achieved 
by requiring that solutions be invariant under the one- dimensional subgroup of the symmetry 
group, corresponding to the Lie algebra element (2.4). In order to decouple these equations 
we had to impose either ai = (no fundamental harmonic), or the supplementary condition 
(3.1) on the phases (phase locked solutions). The existence of the symmetry W in eq.(2.3) 
guarantees that if the phases locking condition (3.1) is imposed on the initial conditions, it 
will survive for all times. 

Thus, by construction, all the explicit solutions obtained in this article are phase locked 
travelling waves. Let us discuss some of their features. 
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Solution (4.14), for certain values of the parameters r, a, f3, 5, characterizing the material 
involved, and of the constants u, v, lo\ and Ki, characterizing initial conditions, can be 
periodic finite waves (see eq.(4.5)). The second harmonic R2 then oscillates in the interval 
( — W^f 1 ' ~ I§^ 2 ) anc ^ the fundamental wave also oscillates between finite limits. 

For other conditions (see eq.(4.18)) the elliptic function solutions reduce to solitary waves 
(4.19) with the second harmonic going through a zero when the fundamental one reaches its 
maximum value. Asymptotically both waves (4.19) tend to zero with the second harmonic 
vanishing at a faster rate. 

Many of the obtained solutions are singular, either at some specific point £ = £0 as in 
eq.(4.20), or periodically, as in (4.7), or (4.11). 

For conditions leading to eq.(3.25) the second harmonic differs from the fundamental 
one just by a proportionality factor and a constant shift of the amplitude (see eq.(3.26)). 
Similarly, for conditions (3.28), the two waves are simply proportional. 

Singular solutions coexist with the finite ones for all the values of the coefficients in 
eq.(2.1). Their physical meaning needs a separate investigation: they may be an indication 
that higher harmonics or dissipative effects were unjustifiably ignored in the derivation of 
these equations. These effects would tend to smoothen out the singularities and possibly 
turn them into finite resonance phenomena. 

The existence of families of elliptic function solutions can be viewed as a manifestation 
of "partial integrability" . We have shown that the studied equations do not, in general have 
the Painleve property. For special values of the constants involved we do get solutions that 
do have this property: they have no movable singularities other than poles. 

Finally let us mention that for particular values of the constants in eq.(2.1) the sym- 
metry group may be larger. For instance, for (3 = 5 = the equations are invariant under 
dilatations generated by: 

D = xd x + 2td t - 2(R 1 d Rl + R 2 d R2 ). 
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This raises the posibility of obtaining 



self-similar solutions of the form 



«i(0 




1 



F 2 (Oe iMO e i 



t 




Since self-similar solutions are particularly important (and stable) in optical systems with 
memory [20,21,22], this situation may be will worth exploring. 
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